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The energy dependence of the total hadronic cross section at high energies is investigated with 
focus on the recent experimental result by the TOTEM Collaboration at 7 TeV and the Froissart- 
Martin bound. On the basis of a class of analytical parametrization with the exponent 7 in the 
leading logarithm contribution as a free parameter, different variants of fits to pp and pp total cross 
section data above 5 GeV are developed. Two ensembles are considered, the first comprising data 
up to 1.8 TeV, the second also including the data collected at 7 TeV. We shown that in all fit 
variants applied to the first ensemble the exponent is statistically consistent with 7 = 2. Applied 
to the second ensemble, however, the same variants yield 7's above 2, a result already obtained in 
two other analysis, by U. Amaldi et al. and by the UA4/2 Collaboration. As recently discussed by 
Ya. I. Azimov, this faster-than-squared-logarithm rise does not necessarily violate unitarity. Our 
results suggest that the energy dependence of the hadronic total cross section at high energies still 
constitute an open problem. 
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I. INTRODUCTION 
A. General Aspects 

High-energy particle collisions constitute the main experimental tool in the investigation of the 
inner structure of matter. In Particle Physics, high-energy usually means center-of-mass energies 
above 10 rup ~ 10 GeV, where nip is the proton mass^ jlj. Presently, for particle-particle and 
antiparticle-particle collisions, the highest energies reached in accelerators concern proton-proton (pp) 
and antiproton-proton (pp) interactions, corresponding to 7 TeV and ~ 2 TeV, respectively. These 
hadronic processes are expected to be described by the Quantum Chromodynamics (QCD), the non- 
Abelian gauge field theory of the strong interactions [2] . 

As a non-Abelian theory, the gluons (the field quanta) themselves carry a color charge and can 
therefore interact with other gluons^. The dynamical consequence is a running coupling constant 
as', the color charge is small at short distances and large at large distances, leading to two different 
regimes named asymptotic freedom and confinement, respectively [H-Q. In hadron-hadron collisions 
these regimes correspond to two sectors that have been known as hard scattering (small distances and 
large values of the momentum transfer) and soft scattering (large distances and small values of the 
momentum transfer). The confinement barrier, where as ~ 1, is typical of distances of the order of 1 
fm and therefore "peripheral" hadronic collisions correspond to the soft sector. 

The great triumph of QCD concerns the perturbative techniques, successfully applied in the hard 
sector [j]. By contrast, soft interactions, characterized by small values of the momentum transfer, can 
not be treated through these techniques due to the rise of the coupling constant as the momentum 
transfer decreases. As a consequence, soft physics demands first principles and nonperturbative ap- 
proaches, which means the non-trivial investigation of the vacuum structure, intricate Monte Carlo 
simulations and complex analytic formalisms However, despite the success of nonperturbative 
QCD in the investigation of the static hadronic properties (bound states), a formal approach to soft 
scattering states, based on first principles and without model assumptions, is still missing 0, 0] and 
that implies in some fundamental problems. 

Soft scattering embodies elastic collisions and diffraction dissociation (single and double) @ and here 
comes one of the striking features of QCD: elastic scattering, the simplest kinematic collision process, 
just constitute one of the greatest dynamic problems for the theory of the strong interactions. In this 
respect, the Optical Theorem Q plays a crucial role since it connects the forward elastic scattering 
amplitude with the most important physical quantity characterizing a collision process, namely the 
total cross section Therefore, the lack of a pure QCD result for the elastic amplitude puts serious 

limits in the theoretical investigation of the total hadronic cross section. On the other hand, and more 
important for our purposes, experimental information on the behavior of the total cross section may, 
in principle, be used as input providing new insights in the development of the theory in the soft 
sector (the inverse problem), at least in what concerns the forward elastic scattering amplitude. 

Operationally the total cross section is defined by 

_ Nel + N.nel 

<^t0t — 7 

where C is the luminosity (fiux per unit area) and N^i, Ninei are the rate of elastic and inelastic 
interactions, respectively (scattered fluxes). From the definition, two interpretations emerges for 
cTtot, one statistical, as a probability of interaction (ratio between incident and scattered particles) 
and another geometrical, as an effective area of interaction, usually measured in mb for hadronic 
scattering. 

Experiments indicate that at lower energies, below ~ 2 GeV, atot is characterized by narrow peaks, 
caused by the formation of resonances (bound states). As the energy increases, reaching the scattering 
region, atot decreases very slowly up to ^ 20 GeV and then starts to grow smoothly and monotonically 
(without bumps or dips) , as illustrated by the experimental points in Figure [1] in the case of pp and 



^ As usual in high-energy physics we adopt the system of units in which c = h = 1. Typical units of energy are 1 GeV 

= 10^ eV and 1 TeV = lO^^ eV. 
^ This property contrasts with QED since photons do not have electric charge. 
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pp scattering (see also Figure 41.10 in [10]). The energy dependence of atot is therefore a crucial point 
since, above the resonance region, it is directly related to the soft sector (optical theorem) giving 
information on the dynamics of the elastic interaction 0, 0] . 

B. Purposes and Goals of the Paper 

Although the rise of the total hadron-hadron cross sections at high energies is an experimental fact, 
the exact energy dependence involved has been a long-standing problem. Several phenomenological 
models, with distinct physical pictures and good descriptions of the available data, make different 
asymptotic predictions [a, [HI, and the only general, formal, widely accepted result is the well- 
known Froissart-Martin (FM) bound [H, Q 

crtot{s) < c In^ — , (1) 
So 

where s is the squared center-of-mass energy, sq a constant, and a bound on the first term on the 
right-hand side, 

c < ir/ml w 60 mb, (2) 

has been obtained by Lukaszuk and Martin [l^. Even recently, the foundations of this key result 
in soft hadronic physics have been discussed in the literature |16l - ll8l |. In particular, Azimov has 
presented a short critical review on the assumptions supporting the derivation of the bound LlSj . 

Even depending on the unknown squared mass scale sq, the numerical values associated with the 
FM bound lie far above the existing data for the total cross section. For example, for sq in the 
interval 1-50 GeV^, the bound in ((T][5]) at 1 TeV is of order of 10 b, much larger than the 100 
mb cross sections typically found in experiments at the highest energies. However, the bound also 
comprises a maximum rate of rise for the cross section with the energy, namely the squared logarithm 
behavior at the asymptotic energy region and that is the point we are interested in here. Different 
dependencies have been extensively tested and discussed in the literature, specially in the context of 
amplitude analysis, which are characterized by analytical parametrization for the total cross section 
and fits to forward data. The most common functional forms employed in the highest-energy region 
include different combinations of constant, linear/quadratic logarithm dependencies and s power laws, 
as discussed, for example, i n ilft-i28il and references therein. Nonetheless, like the case of specific 
phenomenological models [1, llll Il2|. these analyses present good descriptions of the available data 
with different functional forms and hence offer different physical pictures. 

New results from the CERN Large Hadron Collider (LHC) and new estimates from the Pierre Auger 
Observatory for the proton-proton total cross sections are expected to shed light on the subject, not 
only by selecting phenomenological models/pictures, but also by providing information on the degree 
of possible saturation of the bound in terms of the energy dependence of the cross section. In fact, at 7 
TeV the first result for the total cross section by the TOTEM'^ Collaboration, a luminosity-dependent 
measurement (29l |. indicates consistency with a In^ s dependence, as predicted ten years ago by the 
parametrization that was ranked highest by the COMPETE^ Collaboration [13, [Hi, also quoted in 
the Review of Particle Physics by the Particle Data Group [l3]- Therefore, these results favor the 
maximum increase rate allowed by the bound. 

By contrast, at least two previous almost-model-independent analyses, based on data ensembles at 
lower energies, have indicated that the exponent in the logarithm may be somewhat larger than 2 
[30I [3l| . Moreover, on the theoretical side, as recently discussed by Azimov [l^, it is not obvious if 
Martin's derivation, in the context of axiomatic local quantum field theory, can be directly applied to 
hadronic processes. Formal arguments suggest that the total cross section could grow faster than In^ s. 
A recent result from lattice QCD [s^] indicates an universal asymptotic squared logarithm dependence 
for the hadronic total cross section. That conclusion nonetheless rests on specific assumptions and is 
hence neither unique or exclusive. 



^ TOTal Elastic and diffractive cross section Measurement 

* Computerized Models, Parameter Evaluation for Theory and Experiment 
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In view of these facts the fohowing question arises. Taking into account the TOTEM resuh and 
considering the exponent in the logarithm as a free fit parameter, the In^ s dependence is in fact a 
unique solution describing the asymptotic rise of the total cross section, or the data can be statistically 
described by another solution, rising faster (or slower) than In^ s? To answer this question is the goal 
of this work. 

For this purpose, we shall revisit the analytical parametrization introduced by Amaldi et al. in 
the seventies [30|, also employed by the UA4/2 Collaboration, in the nineties [3l|, characterized 
by Reggeon contributions at low energies and a leading contribution at high energies parametrized 
by a power law in Ins, the real exponent 7 as a free parameter. As will be recalled along the 
text, the two analyses cover different energy intervals and led to an exponent exceeding 2. For 
future reference, we note that the above-mentioned highest-ranking parametrization, obtained by the 
COMPETE Collaboration [ll, HE] , can be regarded as a particular instance of this parametrization 
in which the exponent was fixed at 7 = 2. 

All these works developed simultaneous fits to atot and the p parameter (ratio between the real 
and imaginary parts of the forward amplitude) . Here, for reasons to be discussed in detail, we shall 
analize only the Utot data, from pp and pp scattering above -y/s = 5 GeV. In order to investigate the 
effect of the recent experimental result for the pp total cross sections by the TOTEM Collaboration, 
we first consider all the previous existing data, covering the region up to 1.8 TeV, and then add to 
the ensemble the data at 7 TeV. We choose two different initial values for the free fit parameters and 
six fitting variants. We show that, with data up to 1.8 TeV, the real exponent in the logarithm term 
is statistically consistent with 7 = 2, as predicted by the COMPETE Collaboration. However, with 
the addition of the data at 7 TeV, the fits indicate exponents larger than 2 in all cases investigated. 

The manuscript is organized as follows. In Sect. |ll]we introduce the analytical parametrization, 
followed by some comments on its applicability and a critical discussion on the selected data ensemble. 
In Sect. Illll we present the fit procedures, variants and results. The conclusions and some final remarks 
are the contents of Sect. |TVl 



II. ANALYTICAL PARAMETRIZATION AND DATA ENSEMBLE 



A. Analytical Parametrization and Previous Results 

The class of analytical parametrization [s^, 'sT'l consists of two components for the total hadronic 
cross section associated with low- {L) and high-energy {H) contributions: 

CTtot (s) = CTLE {s) + (J HE (s) • (3) 

The first term accounts for the decreasing of the total cross section and the differences between 
particle-particle and particle-antiparticle scattering at low energies and is expressed by 



a-LE{s) 







-bi 






s 


s 


ai 




+ Ta2 













-b2 



(4) 



where si — 1 GeV^ (fixed) while ai, bi, 02, &2 are free fit parameters, and 

-1 for particle-particle 
-1 for antiparticle-particle. 

The second term accounts for the rising of the cross section at higher energies and is given by 

aHE{s)^a +13 In'' (5) 

Sh 

where a, l3, 7, Sh are real free parameters. For further discussion, let us briefly recall some formal 
aspects and previous results associated with this class of parametrization. 

In the context of the Regge-Gribov theory all terms in Eqs. (|4|) and ^ have specific physical 
interpretations, namely Reggeon and Pomeron exchanges at low and high energy regions, respectively 
[25| . The Reggeons correspond to mesons resonances families with the adequate quantum numbers 
in the i-channel process and represented by trajectories interpolating the data on plots of spin J 
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versus the square of their masses (Chew-Frautschi plot). In this case, bi and 62 correspond to the 
intercept of the trajectories and ai, 02 the Reggeon strengths (residues). For pp and pp scattering, 
aLE is associated with two Reggeons, the first one with C = + 1 (02 and /2 mesons trajectories) 
and the second with C = - 1 (p and w mesons trajectories), corresponding to t = + 1 and r — -1, 
respectively. The type of Pomeron contribution depends on the 7 value. For 7 = 1, the constant plus 
Ins terms correspond to a double pole at J = 1 and for 7 = 2 a triple pole (expressed by In^ s. Ins 
and the constant terms). Up to our knowledge, the case of real exponent and < 7 < 2 corresponds 
to a strong-coupling scenario (critical Pomeron) ^33', '34] and a fractional power, 7 = 3/2 (in general 
1 < 7 < 2), is indicated by the mini-jet QCD model with infrared gluon resummation [35l. |3(||. 

As commented before, the above parametrization has been introduced by Amaldi et al. ^30], with 
ule expressed as function of the lab energy, a he as function of s and with a fixed scale constant Sh 
= 1 GeV^ . Simultaneous fits to utot and p data (via dispersion relations) from pp and pp scattering 
in the interval 5 < ^/s < 62 GeV has lead to the result 

7 = 2.10 ± 0.10. 

The same functional form was subsequently used by the UA4/2 Collaboration 'si']. Simultaneous 
fits to Utot and p data from pp and pp scattering in a larger interval, 5 < y/s < 546 GeV, have yielded 
the result 

7 — Z.ZO^Q 

More recently the COMPETE Collaboration has carried out a detailed and extensive study on 
possible analytic parametrizations including all measured Utot and p data from pp, pp, meson-p scat- 
tering, among other processes, at yfs > 4 GeV and pp data up to 1.8 TeV [l^, [l^- Different aspects 
of fit qualities have been considered in a ranking procedure with the same ule structure and a he 
parametrized either with 7 = 1 or 7 = 2 and the constant term. The parametrization with 7 = 2 was 
ranked highest. As commented before, the predictions from this analysis, carried out 10 years ago, 
indicate consistency with the TOTEM result at 7 TeV and therefore agrees with the saturation of the 
squared logarithm dependence in the FM bound. 



B. Data Ensemble and Critical Comments 



Two aspects of our choice of data ensemble deserve special atention. One aspect concerns the 
reactions and the other the physical quantities to be investigated, as discussed in what follows. 

First, the squared logarithm dependence in the FM bound is an asymptotic result, and the 
parametrization ([3][5]) covers particle-particle and particle-antiparticle interactions. For those rea- 
sons we shall consider only the cases with the highest energy interval in terms of available data, 
namely pp and pp scattering. With this restrictive choice we do not take account of any constraint 
dictated by a supposed universal behavior, or data from other reactions in the region of intermediated 
and low energies, as for example the meson-proton cases. Following Amaldi et al. and the UA4/2 
Collaboration, we focus our analysis on data at > 5 GeV^. 

Second, as we have commented before, amplitude analyses of the growth of the total cross section 
include information on the p parameter, through dispersion relations (integral and/or derivatives 
forms), or the asymptotic prescriptions for crossing even and odd amplitudes (Phragmen-Lindeloff 
theorem) (37. .39] . With the exponent 7 as a real free fitting parameter, the integral dispersion relations 
demand numerical methods and therefore require an specific approach for error propagation from the 
uncertainty in the 7 parameter. The use of prescriptions seem to us unjustified in the region of 
intermediate and low energy data, since they are asymptotic results (STj . By constrast, derivative 
dispersion relations allow an analytical approach and can be extended down to 4 - 5 GeV [13, Slj or 
even below (above the physical threshold) in the form of a double infinite series [1^, or a single 
series [IJ. However, and that is a crucial point in this work, we shall not consider simultaneus fits to 
total cross section data and p information for the six reasons that follows. 

1. Strictly speaking the p parameter is not a quantity with the same physical status as the total 
cross section since, in practice, it is evaluated as a free fit parameter in the Coulomb-nuclear 
interference region or inferred from analytical parametrization [38l . 139^ . 

2. As the energy grows, it becomes progressively more difficult to determine p [H,!!^ and therefore 
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the associated uncertainty, as can be easily seen in plots of this parameter in terms of the energy. 
Even the COMPETE Collaboration refers to the difficulty to adequately fit the p data from pp 



3. Simultaneous fits to total cross section data and p parameter demand the use of dispersion rela- 
tions with one subtraction psl l46l | and therefore one more parameter, the subtraction constant. 
However this constant does not have a physical interpretation, constrasting therefore with all 
the parameters in the atot parametrization (shortly discussed in Section [II A[) . 

4. In data reductions, the correlation of the subtraction constant with all the other (physical) 
parameters affects the fit results at both low- and high-energy regions, as discussed in detail in 
[43] ■ As a consequence, the presence or not of the subtraction constant may lead to different 
results. In this respect we note that, although referred to in [sO] the value of the subtraction 
constant is not given (null or neglected?), in [31j its fit value is — 57 ± 4 and it is neglected in the 
prescriptions or derivative dispersion relations used by the COMPETE Collaboration [U,!!^. It 
should also be noted that even the prediction of the p parameter from fits to total cross section 
data is affected by the presence or not of the subtraction constant (47l . fisj . 

5. As recently demonstrated by Ferreira [i^ and collaborators [13, [Ell (see also ^55]), the different 
slopes associated with the real and imaginary parts of the hadronic amplitude in the Coulomb- 
nuclear interference region, affect the extracted value of the p parameter. This effect, however, 
has never been considered in the experimental procedures used in the p determination. The 
limited validity of the relative phase between Coulomb and nuclear amplitudes, used in the 
experiinental procedures to determine p, has also been discussed by Kundrat, Lokaji'cek and 
Vrkoc [53]. 

6. As a consequence of the above mentioned effects and the increasing uncertainty in the p de- 
termination with the energy, any possible deviation from an analytical parametrization for the 
total cross section, dictated by the experimental data at the highest energies, may be hidden or 
lost. In other words, the inclusion of the p information may erroneously anchor the rise of the 
total cross section at the highest energies. 

These six critical points suggest that, at the high-energy region, fits restricted to the total cross 
section data may avoid the bias introduced by both the p parameter and the subtraction constant 
embodied in the dispersion relations. Therefore, although not being a usual procedure in amplitude 
analysis, we understand that to explore the possibility of focusing our fits only on the total cross 
section data constitute, at least, a valid strategy. For that reason our data sample comprises only the 
pp and pp total cross section data above 5 GeV including in special, the recent TOTEM result 
at 7 TeV [H. 

Two points must be stressed in the analysis that follows: [i] we are concerned only with the rate of 
increase of the total cross section, not with the numerical value of the EM bound; [ii] our data ensemble 
is only a subset of the data employed in the global analysis by the COMPETE Collaboration. 



By means of a luminosity-dependent method the TOTEM Collaboration has recently obtained for 
the pp total cross section at 7 TeV the value crf^j = 98.3 ± 0.2^'''* ± 2.8''^''^ mb Our point here 
is to investigate the effect of this new data in the previous analytical fits covering the energy interval 
5 GeV < ^/s < 1.8 TeV, as was the case with the COMPETE Collaboration. Based on the results by 
Amaldi et al. and by the UA4/2 Collaboration, our strategy amounts to using the parametrization 
dSKS]) to investigate possible deviations from 7 = 2. 

However, several aspects related to fitting procedure and the alternatives allowed by physical con- 
siderations or assumptions deserve investigation. The three main points concern: (1) the nonlinearity 
of the fit demands a methodology for the choice of the initial (feedback) values of the free parameters 
(53 |: (2) as commented in Sect. Ill A[ in the theoretical context the intercepts hi and ^2 are expected to 
be consistent with spectroscopic data (Chew-Frautschi plot); (3) fixing hi and 62 affects the fit results 
in both low- and high-energy regions, due to the correlation among the free parameters. In order to 
address these points the following methodology and fit variants have been considered. 




III. FITTTING PROCEDURE, VARIANTS AND RESULTS 
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In all cases under study we first consider the ensemble in the interval 5 GeV < ^/s < 1.8 TeV and 
after that we include the data at 7 TeV from the TOTEM Collaboration. For easier reference we 
will denote these two ensembles by \/s^^^ — 1.8 TeV and -s/s^a^. = 7 TeV, respectively. In order to 
investigate the effects of different choices for the initial values of the free parameters, we have followed 
two alternative procedures, named Method 1 and Method 2: 

Method 1. 

Since the highest rank result by the COMPETE Collaboration predicts a In^ s dependence and is 
consistent with the TOTEM result at 7 TeV, we consider as a kind of orthodox choice to initialize our 
parametric set with the central values they have obtained in the simultaneous fit to atot and p data, 
which includes pp, pp in the interval 5 GeV < i/s < 1.8 TeV as a subset. The numerical values [l3| 
are displayed in the second column of Table HI 

Method 2. 

Alternatively we have chosen bi — b2 ~ 0.5 (average values for Reggeon intercepts), 7 = 2, ai = 02 
= a = 1 mb, s/i = 1 GeV^ and (3 — 50 mb (simulating a saturation of the Lukaszuk-Martin bound 
[20I ) . The numerical values are displayed in the first column of Table IIIII 

The data reductions were carried out with the objects of the class TMinuit of the ROOT Framework 
[ssl ]. A Confidence Level (CL) of w 68 % (one standard deviation), was adopted in all fits, so that the 
projection of the distribution in (A^ + l)-dimensional space (A^ = number of free fit parameters) 
has the probability of 68 % ^54] . 

With both methods and the two ensembles, different variants were also considered. Sections IIIIAI 
and nil B] describe the variants and present the results obtained with Methods 1 and 2, respectively, 
the discussion of the results being deferred to Sect. |TVl 

A. Method 1 

1. Results for the \/Smai ~ 1-8 TeV Ensemble 

We consider the cases and notation that follows. The numerical results and statistical information 
are displayed in Table HI 

Direct Fit: Initialized with the COMPETE parameters and 7 = 2 (fixed). The first run of the 
MINUIT Code yields the for that ensemble (second column in Table H]) and the final run 
gives the direct fit result for that ensemble (third column in Table |I| . 

Variant 1 (VI): Initialized with the Direct Fit parameters but now with the exponent 7 as a free 
parameter (fourth column in Table IJ). 

Variant 2 (V2): Also initialized with the Direct Fit results. In this case, we consider 7 = 2 (fixed) 
and hi and 62 also fixed to the intercepts extracted from the spectroscopic data on the 02 / /2 
and p/lu mesons trajectories, obtained by Luna, Menon and Montanha, namely hi = 0.452 and 
&2 0.558 [56^ (fifth column in Table U). 

Variant 3 (V3): Also initialized with the Direct Fit results. In this case we consider 7 free and 61 
and 62 fixed as in Variant 2 (sixth column in Table |l|. 

In terms of the initial values and data reductions the following diagram summarizes the fitting 
procedures: 



{Variant 1 
Variant 2 
Variant 3 
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2. Results for the y/s^^^ = 7 TeV Ensemble 

Here we consider as initial values each one of the corresponding results obtained with the \/s^^^ 
= 1.8 TeV ensemble, namely those listed in Table HI The numerical results and statistical information 
are displayed in Table |TT1 

TABLE I: Results of Method 1 for the \/s„^^ = 1.8 TeV ensemble. Fit result from the COMPETE parameters 
as initial values, Direct Fit and Variants 1, 2 and 3 (initial values from the Direct Fit result), together with 
the statistical information: degrees of freedom (DOF) and reduced x^- The parameters ai, a2, a and /3 are in 
mb, Sh is in GeV^ and 6i, &2 and 7 are dimensionless. 





Initial Values 


Direct Fit 


VI 


V2 


V3 




(COMPETE) 










ffll 


42.53±0.23 


54.6±4.0 


54.6±1.7 


55.4±3.1 


57.0±2.5 


bi 


0.458±0.017 


0.491±0.067 0.4907±0.0095 0.452 (fixed) 


0.452 (fixed) 


a2 


33.34±0.033 


33.1±2.3 


33.1±1.7 


35.78±0.36 


35.78±0.39 


62 


0.545±0.007 


0.540±0.016 


0.540±0.012 


0.558 (fixed) 


0.558 (fixed) 


a 


35.45±0.48 


34.2±2.5 


34.19±0.24 


32.14±0.99 


31.38±0.98 


P 


0.308±0.010 


0.264±0.029 


0.263±0.018 


0.245±0.018 


0.180±0.062 


7 


2 


2 (fixed) 


2.001±0.026 


2 (fixed) 


2.083±0.095 


Sh 


28.9±5.4 


12±12 


12.2±1.5 


5.7±3.0 


2.8±2.6 


DOF 


156 


156 


155 


158 


157 


xVdof 


1.02 


0.931 


0.937 


0.929 


0.934 



TABLE IL Results of Method 1 for the y/s^^^ = 7 TeV ensemble. Fit results with initial values corresponding 
to the final values given in Table HI Direct Fit and Variants 1, 2 and 3. Units as in Table H) 





Direct Fit 


VI 


V2 


V3 


ai 


54.9±5.3 


57.7±2.0 


53.8±3.0 


58.4±2.4 


bi 


0.539±0.080 


0.526±0.010 


0.452 (fixed) 


0.452 (fixed) 


02 


33.2±2.3 


33.2±1.7 


35.77±0.36 


35.78±0.36 


b2 


0.541±0.016 


0.541±0.012 


0.558 (fixed) 


0.558 (fixed) 


a 


35.9±1.9 


34.90±0.23 


32.68±0.87 


30.41±1.00 


P 


0.290±0.027 


0.199±0.014 


0.257±0.017 


0.093±0.019 


7 


2 (fixed) 


2.104±0.027 


2 (fixed) 


2.273±0.039 


Sh 


26±20 


11. Oil. 4 


7.8±3.5 


0.77±0.58 


DOF 


157 


156 


159 


158 


XVI50F 


0.930 


0.935 


0.931 


0.933 



The corresponding curves obtained with Method 1 for both ensembles and the experimental data 
are shown in Figures [T] to 21 In each figure the COMPETE result is displayed as reference (solid lines), 
together with each result obtained with ensembles = 1-8 TeV (dot-dashed lines) and 

— 7 TeV (dotted lines): Direct Fit and Variants 1, 2 and 3 in Figures [l] [2l [3] and HI respectively. 

In these figures it is also displayed estimates for the total cross section from cosmic-ray experiments 
at the highest energies, which were not included in the data reductions. These estimates, discussed 
in [l^l (see also were obtained by the Akeno Collaboration [5^ in the ~ 6 - 24 TeV region, and 

by the Fly's Eye Collaboration at 30 TeV [58]. Also included is the recent result by the Pierre Auger 
Collaboration at 57 TeV 59]. We shall discuss these and the results that follows in Section HVl 
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B. Method 2 



1. Results for the \/Smai^ ~ ^-^ TeV Ensemble 

The initial values, referred to in the beginning of this Section, are displayed in the second column 
of Table Hill Here, the following variants have been considered: 

Variant 4 (V4): The parameters 6i and 62 have been fixed to the average values expected for the 
Reggeon intercepts and 7 = 2 also fixed (third column in Table HlH) . 

Variant 5 (V5): Same as Variant 4 with 7 as free parameter (fourth column in Table Hill . 

Variant 6 (V6): Initialized with the parameters from Variant 4: 61, &2 and 7 free parameters (fifth 
column in Table HlH) . 

The following diagram summarizes the fit procedures and initial values (Table IIIip : 



Initial Values 



Variant 4 — > Variant 6 
Variant 5 



2. Results for the y/Srnax = T TeV Ensemble 

Here we followed the same procedure described above, now with the expanded ensemble. Table HVl 
displays the numerical results and statistical information. Figures [5l [5] and [7] compare the accelerator 
data and estimates from cosmic-ray experiments with the fits obtained with Method 2 and variants 4, 
5 and 6, respectively for the \/s^^^^^ = 1.8 TeV (dot-dashed lines) and \/s,„aa: — ^eV (dotted lines) 
ensembles. 

TABLE III: Results of Method 2 for the = 1.8 TeV ensemble. Fit results with Variants 4, 5 and 6. 

Units as in Table HI 





Initial Values 


V4 


V5 


V6 


ai 


1 


52.6±3.4 


54.7±2.6 


53.8±1.8 


61 


0.5 


0.5 (fixed) 


0.5 (fixed) 


0.4952±0.0099 


0,2 


1 


27.70±0.28 


27.71±0.28 


33.1±1.7 


62 


0.5 


0.5 (fixed) 


0.5 (fixed) 


0.540±0.012 


a 


1 


34.86±0.69 


34.32±0.77 


34.60±0.24 


P 


50 


0.270±0.019 


0.21±0.11 


0.290±0.021 


7 


2 


2 (fixed) 


2.06±0.16 


1.975±0.028 


Sh 


1 


15.6±6.3 


9.3±8.4 


16.0±2.0 


DOF 




158 


157 


155 


xVdof 




0.963 


0.969 


0.937 



We note a difference in the diagrams related to initial values and variants (two schemes displayed in 
this Section) . This effect is consequence of fit procedures since in certain cases the data reductions led 
back to the initial values or to solutions with errors that exceded the central values of the parameters. 
All these cases have been excluded. 



IV. CONCLUSIONS AND FINAL REMARKS 



We have employed two methods to initialize the fitting parameters, six variants and two data 
ensembles. The results for the \fs^ax = 1-8 TeV ensemble, listed in Tables HI and IIIIl are statistically 
consistent with 7 = 2, which points to a saturation of the squared logarithm dependence for the total 
cross section at high energies. More specifically, the results can be summarized as follows: 
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TABLE IV: Results of Method 2 for the a/s„„^ = 7 TeV ensemble. Fit results with Variants 4, 5 and 6. Same 
units as in Table H] 





V4 


V5 


V6 


ai 


51.34±3.2 


56.5±1.1 


57.6±7.5 


bi 


0.5 (fixed) 


0.5 (fixed) 


0.525±0.077 


a2 


27.70±0.28 


27.70±0.28 


33.2±2.3 


b2 


0.5 (fixed) 


0.5 (fixed) 


0.541±0.016 


a 


35.13±0.60 


33.65±0.22 


34.9±1.7 


P 


0.279±0.016 


0.1301±0.0086 


0.199±0.064 


7 


2 (fixed) 


2.213±0.024 


2.10±0.11 


Sh 


18.5±6.3 


3.90±0.52 


10.8±5.6 


DOF 


159 


158 


156 


xVdof 


0.962 


0.967 


0.935 



Method 1 (Table U: 7 « 2.00 ± 0.03 (VI) and 7 « 2.08 ± 0.10 (V3). 
Method 2 (Table ini|: 7 « 2.06 ± 0.16 (V5) and 7 w 1.98 ± 0.03 (V6). 

On the other hand, when the more recent 7 TeV TOTEM data is included, all fits with 7 as a free 
parameter, listed in Tables HIl and IIVI are statistically consistent with exponents above 2, as indicated 
in the following summary: 

Method 1 (Table HI): 7 ~ 2.10 ± 0.03 (VI) and 7 w 2.27 ± 0.04 (V3). 
Method 2 (Table HVll: 7 « 2.21 ± 0.02 (V5) and 7 « 2.10 ± 0.11 (V6). 

All these results for 7, together with those obtained by Amaldi et al. and the UA/4 Collaboration, 
are schematically displayed in Figure [S] From our data reductions through parametrization ([3][5]) to 
pp and pp scattering above 5 GeV, including the 7 TeV TOTEM result, we conclude that the total 
hadronic cross section may rise faster than In^ s at high energies. 

From Figure [HI we see that the central values of the 7 parameter from the fits with \/s^^^^ = 1.8 
TeV (Methods 1 and 2) are below the corresponding values obtained by Amaldi et al. and by the 
UA4/2 Collaboration. Since these authors have included the p information in their analysis, it seems 
not clear that simultaneous fits to p and atot data may, in that case, anchor the rise of the total cross 
section. In other words, we understand that the results we have obtained with v/Smaa: = 1-8 TeV and 
V^maa: ~ TcV cnsembles are not connected with the inclusion or not of the p information, but are 
direct consequences of the total cross section data analyzed and the possibility to treat 7 as a free 
fit parameter (not fixed to 2). In fact, our predictions for p{s) in the case of 7 > 2 are presented in 
Appendix |2] and show good agreement with experimental information available above 5 GeV. 

Concerning the rise of the total cross section faster than the squared logarithm dependence, we 
have the following final comments: 

1. In the theoretical context, different assumptions on the nonphysical amplitudes in the asymptotic 
region can explain this behavior, which does not mean violation of unitarity, as recently discussed by 
Azimov [13. 

2. In the experimental context, beyond the LHC energy region, extremely large uncertainties 
are associated with estimations of the pp total cross section from cosmic-ray experiments. Such 
uncertainties are unavoidable since the flux decreases as the energy grows and because the extraction of 
(Tpp from <Jp-air is model dependent. These estimations have been included in our figures for qualitative 



comparison only. It may be useful to recall that the the results of the Akeno Collaboration [57[ have 



been criticized in several works, as discussed for example in Refs. [25:| and 47[, and references therein. 
Had we ignored this information, the Fly's Eye and Auger results, i.e., the highest-energy points in our 
figures, might suggest a different scenario for the rise of the total cross section. In fact, from Figures [T] 
to [71 all the curves in consistency with the TOTEM data lie above the central values of these cosmic- 
ray estimations and the same is true in the inverse sense. That is, there is no agreement among the 
TOTEM result and the Fly's Eye and Auger estimations, at least in what concerns parametrization 



D.A. Fagundes, M.J. Menon, P.V.R.G. Silva 



11 



(I3][5]). In this sense, the expected measurement of atot by the TOTEM CoUaboration at 7 and 8 
TeV through the luminosity-independent 0, method may shed hght on the subject. Moreover, new 
experimental results on elastic and diffractive scattering at 8 TeV (and, subsequently, at 14 TeV), 
will provide novel phenomenological insights and reduce the uncertainties from model extrapolations 
necessary to obtain atot from cosmic-ray experiments (see "go"] and references therein). 

At last we stress that our approach and strategies do not follow the usual or standard lines of the 
amplitude analyses on the energy dependence of the total cross section. As commented, the latter 
include p data, associate the parametrization with data at low energies and constrain the fits to 
describe other reactions, for which low-energy data are available. These are aspects that we expect to 
consider in a future work, since they may provide information that is complementary to the results 
we have presented. As commented along the manuscript, our goal has been to explore the possibility 
of investigating only the reactions with highest energy interval available, concentrating on the total 
cross section data. We have tried to identify possible high-energy effects that may be unrelated to the 
trends of the lower-energy data (including those from other reactions) or hidden in fits including the p 
information. Our intention is not to compete with other authors or analysis, but to call the attention 
to the possibility that the rise of the total hadronic cross section may still constitute an open problem, 
an assertion that contrasts with the view advocated by some authors [281] . 

Note added during revision 

After this paper was submitted to publication, new fits of forward quantities were obtained by the 
Particle Data Group The updated data set includes the TOTEM result at 7 TeV and the new 
fit with the highest rank COMPETE parametrization (7 = 2) shows that the point at 7 TeV is not 
described: the curve lies below the data (Figure 46.10 in [6l|). This result corroborates those already 
shown with the y/s„^ax — ^ TeV ensemble in our Figures [3 [3] and [5] 
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Appendix A: Predictions for the ratio between the real and imaginary parts of the forward 

amplitude 

In this appendix we present predictions for p(s) from fits to pp and pp total cross section data through 
parametrization ([SE]) with 7 > 2. We also discuss the effect and role of the subtraction constant, 
embodied in the associated dispersion relation. In search of extreme effects we shall consider the two 
results obtained with the \^jnax ~ TeV ensemble, which correspond to the highest exponents 7, 
namely 

Method 1 - Variant 3: 7 w 2.27 ± 0.04 (Table HI) 
Method 2 - Variant 5: 7 « 2.21 ± 0.02 (Table |IV| 

In terms of the forward scattering amplitude F{s), the total cross section (optical theorem) and the 
ratio between the real and imaginary parts of the amplitude, at high energies, can be expressed by 

, , lmF{s) . . ReF{s) 

crMs) = Pis) = (Al) 

For crossing even (-I-) and odd (— ) amplitudes, pp and pp scattering are related by 

F±{s) = ^IE±IhL (A2) 
and it follows from Eqs. (IXTj) and (|X2|) that 
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pPP 



■Zis) 



ReF4- 



ReF^ 



(A3) 



and 



1 r Re I Re F^ 



(A4) 



The real and imaginary parts of the even and odd amphtudes are connected by dispersion relations 
and the high-energy domain demands one subtraction psl |46| . Since we are looking for analytical 
results, we shall work here with the derivative dispersion relations in the standard form deduced by 
Bronzan, Kane and Sukhatme (62l |. They are obtained from the integral dispersion relations in the 
high-energy limit: 



Re F+{s) _ K 
s s 



tan 



2 rflns 



Im F+(s) 



Re F_(s) 



= tan 



2 \ ' dins 



Im F_ (s) 



where K is the subtraction constant. Operationally these relations can be evaluated through the 
expansions [gsL IgI] 



Re I£ 
s s 



2 dins 



2 dins 



15 V 2 dins 



-I- 



Im F+(s) 



(A5) 



Re F^ (s) 




Im F^ (s) 



din ; 



(A6) 



With parametrization ([3][5]) taken as input, a closed form results from the sum of the contributions 
associated with the power-law term gle- The sum of the contributions from the logarithm term a he 
converges fast for 7 in the range (2.2 - 2.3) and a third-order approximation is therefore sufhcient. 
We obtain for the even part 



ReF+(s) 



K [nhi 

fli tan 

s I 2 



) 


s 







n -bi 



-3 



Cln'' 



-5 



^In'^-i 



(A7) 



where 



B 



(7[7-l][7-2], 



C = — 



/?7[7-l][7-2][7-3][7-4] 



(A8) 



and for the odd part 
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Rei^_(s) 




) 



r 1 -b2 
s 



(A9) 



s 



With the above resuhs, Eqs. (jA3[) and (IA4[) yield analytical expressions for pPP{s) and pPP{s). Note 
that the analytical results imply that as s — ^ oo: 



In what follows we pre sent the predictions for p(s) and compare the results with the experimental 
information available [Ifll- To study the effect of the subtraction constant, we follow two alternatives: 
in the first we fix it at = 0; in the second we let ivT be a fit parameter to the p data only. More 
specifically, we take as input the parameters in Tables |TT] and llVi with si — 1 GeV^ . In the first 
alternative, K — 0, we obtain the direct prediction for p{s). In the second alternative we fix all the 
other parameters and adjust K to fit the p data. In the latter case, the first run of the MINUIT Code 
yields the x^/DOF for K = 0. The predictions with the results by Method 1 - Variant 3 (Table HU 
are displayed in Figure IH] and those with the results by Method 2 - Variant 5 (Table IIVI) in Figure [TUl 
We are led to the following conclusions: 

1. Even in these extrema cases, with 7 : 2.2 — 2.3, all the experimental information on p^^{s) and 

p^^{s) above 5 GeV are quite well described; 

2. The subtraction constant affects the low-energy results, below y/s ^ 20 GeV; 

3. The best x^/DOF results (closest to 1) are obtained with if as a free fit parameter. Obviously, that 

is a consequence of adding one more free parameter, however without the physical interpretation 
associated with those present in the parametrization of the total cross section. 

Finally, we recall that in simultaneous fit to atot and p the subtraction constant affects both the low- 
and high-energy regions [13,13. That is a consequence of the strong correlation among the subtraction 
constant and all the other physical free fit parameters. We plan to discuss this consequence and other 
aspects of the fit procedures in a forthcoming paper. 
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FIG. 1: (Color online) Comparison between experimental data and fit results. The dot-dashed and dotted 
lines display the results of the Direct Fit-Method 1 applied to the y/s^^^ = 1.8 TeV and \/^maa: ~ ''' TeV 
ensembles, respectively. The solid line is the fit obtained by the COMPETE Collaboration. Numerical 
information displayed in Table U (second and third columns) and Table lU (second column). 
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FIG. 2: (Color online) Analogous to Fig. [T] the dot-dashed and dotted lines computed with Variant 1 of 
Method 1 (VI). Numerical information displyed in Tables |I] (fourth column) and HIl fthird column). 
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FIG. 3: (Color online) Analogous to Fig. \T\ the dot-dashed and dotted lines computed with Variant 2 of 
Method 1 (V2). Numerical information displyed in Tables |I] (fifth column) and HIl (fourth column). 




FIG. 4: (Color online) Analogous to Fig. \T\ the dot-dashed and dotted lines computed with Variant 3 of 
Method 1 (V3). Numerical information displyed in Tables HI and HTl flast columns). 
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FIG. 5: (Color online) Analogous to Fig. [T] the dot-dashed and dotted lines computed with Variant 4 of 
Method 2 (V4). Numerical information displayed in Tables Hill fthird column) and IIVI (second column). 
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FIG. 6: (Color online) Analogous to Fig. [T] the dot-dashed and dotted lines computed with Variant 5 of 
Method 2 (V5). Numerical information displyed in Tables Hill ffourth column) and IIVI (third column). 
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FIG. 7: (Color online) Analogous to Fig. \T\ the dot-dashed and dotted lines computed with Variant 6 of 
Method 2 (V6). Numerical information displayed in Tables Hill ffifth column) and IIVI (fourth column). 
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FIG. 8: Results for the exponent 7 as a free parameter in different data reductions through parametrization 
Pl5)) . Shown are the results by Amaldi et al. ^ (V^max = 62 GeV), the UA4/2 Collaboration HH 
— 546 GeV), and from our analyses for the v^^maa; = 1-8 TeV and v^^^^. = 7 TeV ensembles with Methods 
1 and 2. 
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FIG. 9: (Color online) Comparison between experimental results for p and the analytical expressions derived 
in Appendix [X] with the parameters obtained by Variant 3 of Method 1, listed in Table [III which yield 
7 « 2.27 ± 0.04. The filled and the open circles represent the experimental data for pp and pp scattering 
respectively. The top panel shows our predictions with K = 0, which correspond to / DOF = 1.93. The 
bottom panel, is the best fit to the p data with K a free parameter {x^ /DOF = 1.45). 
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FIG. 10: (Color online) Analogous to Fig. [9l the analytical results plotted for the parameters obtained with 
Variant 5 of Method 2 (Table HV)) . which yield 7 « 2.21 ± 0.02. The ratios x^/DOF in the top and bottom 
panels are 2.70 and 1.58, respectively. 



